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Abstract 

Unlike the Wilson line in QED the Wilson line in QCD contains path ordering. In this paper 
we get rid of the path ordering in the light-like Wilson line in QCD by simplifying all the infinite 
number of non-commuting terms in the SU(3) pure gauge. We prove that the light-like Wilson line 
in QCD naturally emerges when path integral formulation of QCD is used to prove factorization 
of soft and collinear divergences at all order in coupling constant in QCD processes at high energy 
colliders. 
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I. INTRODUCTION 


In Feynman diagrams the infrared divergences appear whenever the energy-momentnm 
involved with the massless particle becomes very small. Similarly the collinear divergences 
occur when the momenta k,p of two massless particles become parallel in the region 0 < 
k « p. Typically the soft and collinear divergences occur in the Feynman diagrams due to 
momentum integration in the quantum loop diagrams involving massless propagators and 
due to momentum integration in the Feynman diagrams involving emission/absorption of 
massless particles. In quantum electrodynamics (QED) the massless particle is photon and 
in quantum chromodynamics (QCD) the massless particle is gluon. The soft and collinear 
divergences are severe in QCD than that in QED because massless gluons interact with each 
other whereas massless photons do not interact with each other. Since massless particle is 
always light-like one hnds that soft and collinear divergences can be described by light-like 
Wilson line. 

However, the physical quantities measured are all soft and collinear divergences free. 
Hence it is important to prove that all the non-canceling soft and collinear divergences 
in the perturbative Feynman diagrams are factorized in the dehnition of the (physical) 
gauge invariant non-perturbative quantities in QCD such as in the dehnition of the parton 
distribution function and fragmentation function at high energy colliders because the soft 
and collinear limit corresponds to long distance regime. This is done by supplying Wilson 
line in the dehnition of the parton distribution function and fragmentation function |l|. The 
factorization refers to separation of short-distance ehects from the long-distance ehects in 
quantum held theory. 

The proof of factorization theorem in QCD is very non-trivial by using the diagrammatic 
method of QCD 2|, yj but it is enormously simplihed by using the path integral method of 
QCD { 4 , The main idea behind the path integral method of QCD to prove factorization 
is to study the soft and collinear behavior of non-perturbative correlation function such 
as < > in QCD due to the presence of light-like Wilson line 

in QCD. Note that a light-like quark with light-like four-velocity produces SU(3) pure 
gauge potential at all the time-space points x^ except at the ^atial position x transverse 
to the motion of the quark at the time of closest approach j^, 6, 7]. The soft and collinear 
divergences in Feynman diagrams in QCD can be studied by using Eikonal approximation 
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for the propagators and vertices 


QQ,!^ 


15j. Hence dne to the Eikonal approximation for 


soft and collinear divergences arising from the soft and collinear glnons interaction with the 
light-like qnark, the light-like qnark finds the glnon held as SU(3) pnre gange [4, l5|. 

The U(l) pnre gange 


A'^{x) = df^oj{x) 


( 1 ) 


gives the light-like Wilson line in QED 


ie dx^Afj,(x) 


( 2 ) 


which is nsed to stndy factorization of soft and collinear divergences in QED [8|, ll3|. In 
QCD the SU(3) pnre gange 


= -[d^U{x)]U-\x), 
W 


U{x) = 


( 3 ) 


gives the light-like Wilson line in QCD 


Ve 


i9T°- /J-'’ dxi^Al(x) 


( 4 ) 


which is used to study factorization of soft and collinear divergences in QCD 1^. Note 
that, unlike the Wilson line in QED in eq. ([2]) which does not contain path ordering P, the 
Wilson line in QCD in eq. (jl]) contains path ordering V. 

In this paper we get rid of the path ordering V in the light-like Wilson line in QCD by 
simplifying all the inhnite number of non-commuting terms in the SU(3) pure gauge in eq. 

We hnd that the light-like Wilson line in QCD without path ordering is given by 


Pexp[i£/ 

J X4 


= exp[igT°-\ 


2l-D[A{xf)] 


dg 


X ex.])[-igT’’[ 


2 I ■ D[A{x,)] 


I ■ 


d[gA{xi)]. 
dg ■ 


( 5 ) 


where the right hand side of the above equation does not contain path ordering V. In eq. 
(15|) the is the covariant derivative, is the light-like four velocity and A^°‘{x) is the 

SU(3) pure range in QCD, which unlike U(l) pure gauge A^[x) in QED, contains inhnite 
powers of p. 
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Since the light¬ 
en the end points 


ike Wilson line in QCD does not depend on the path but depends only 
4 Is] we hnd from eq. ([5]) that the non-abelian phase or the gauge link in 


QCD without path ordering is given by 

p^-i9j;^dXl-A'^ix+lX)T'^ = exp[z^T“[; 


■ 


d[gA{x)]. 


( 6 ) 


'2l-D[A{x)Y dg 

which is used to study factorization of soft and collinear divergences in QCD where the right 
hand side of the above equation does not contain the path ordering V. 

In this paper we will provide a derivation of eq. ([S]). 

In jd] we have shown that the light-like Wilson line in QCD naturally emerges when path 
integral formulation is used to prove NRQCD factorization at all order in coupling constant 


in heavy quarkonium production. Similarly, in we have shown that the light-like Wilson 
line in QCD naturally emerges when path integral formulation is used to prove factorization 
of soft and collinear divergences of the gluon distribution function at high energy colliders 
at all order in coupling constant. In this paper we will prove that the light-like Wilson line 
in QCD naturally emerges when path integral formulation is used to prove factorization 
of soft and collinear divergences of the quark distribution function at high energy colliders 
at all order in coupling constant. Hence we hnd that the light-like Wilson line in QCD 
naturally emerges when path integral formulation of QCD is used to prove factorization of 
soft and collinear divergences at all order in coupling constant in QCD processes at high 
energy colliders. 

The paper is organized as follows. In section II we derive the light-like Wilson line in QCD 
without path ordering as given by eq. (jS]). In section III we study the gauge transformation 
of the light-like Wilson line in QCD without path ordering. In section IV we prove that the 
light-like Wilson line in QCD naturally emerges when path integral formulation of QCD is 
used to prove factorization of soft and collinear divergences at all order in coupling constant 
in QCD processes at high energy colliders. Section V contains conclusions. 


II. LIGHT-LIKE WILSON LINE IN QCD WITHOUT PATH ORDERING 

The SU(3) pure gauge in QCD is given by eq. ([3]) which contains inhnite number of 
non-commuting terms. Simplifying all the inhnite number of non-commuting terms in eq. 
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we find that the SU(3) pure gauge A^^{x) is given by jo] 

' ^gM[x) _ 


= d^uj\x) 


gM{x) 


ab 


where 


( 7 ) 


Mabix) = r’^^^CO^^ix) 


( 8 ) 


Expanding the exponential in eq. ([7]) we hnd 

= \&‘ui\x)] [1 + |m(i) + + ^M\x) + (9) 

In QED the U(l) pure gauge potential produced by a point charge e is linearly proportional 
to the electric charge e 0,y,0, i.e. , 


d^u}{x) oc e. 

Since u{x) is linearly proportional to e we hnd that a;“(x) is linearly proportional to g 
Since u°‘{x) is linearly proportional to g we write 



u^{x) = ( 11 ) 

where ^“(x) is independent of g. Using eq. ffTTl) in ([9]) we hnd 

-hl^“(a;) = [df^/3\x)] [1 + ^N{x) + ^N\x) + ^N^{x) + ...U (12) 

where 


Nat(x) = r‘’‘P‘(x). 


(13) 


Multiplying g‘^Nab{x) in eq. flT^ we obtain 




{g^f 




[gN{x)A^{x)r = [d^(5\x)] [g^N{x) + + ^iV^(x) + ^-^N\x) + 


2 ! 


3! 


4! 


(14) 


Adding d^l3'^{x) in eq. (IT^ we hnd 


2\2 


(gY 


(gY 


2\3 


2\4 


(gY 


D>^[A{x)]n^) = [d^/3\x)] [1 + g^N{x) + ^N\x) + ^N^x) + ^-^N\x) + ...] 


3! 


4! 


ab 


(15) 
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where 




( 16 ) 


Multiplying in eq. 0121) and then taking derivative with respect to we obtain 

= 19-y(x)] 11 + 

Since right hand sides of eqs. ffTSjl and flTTI) are equal we hnd 

DnA(.m^) = 

Converting /3“(x) to a;“(x) by using eq. flTTD we hnd from eq. flTSD 

D>‘lA{x)]u'‘{x) = 

Multiplying the same independent four vector in eq. flT^ we hnd 

I • = 2 1- D[A{x)]uj\x). 

dg 


Dividing I ■ D[A{x)] from left in eq. fl20|) we obtain 


UJ uc = 


'2 / • D[M(x)]' 
which gives the non-abelian phase 


d[l ■ gA\x)] 
dg 


I ■ 


d[gA{x)]. 


$(a;) = = expligT'^l- 


^2 I ■ D[A{x)] dg 

i\9A(x)]„. 


I- 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


^2 I ■ D[A{x)] dg 

From jd, 51 we hnd that the light-like Wilson line in QCD for soft and collinear divergences 


is given by 

dxi^A-{x)T- ^ ^igT-u9-{xf)^-igT>>ujf>{xi) ^ 


/"<iAZ-A“(x/+ZA)r“j pgig d\l-A>>{xi+l\)T>> 


(23) 


Using eq. (|2^ in eq. (|23l) we hnd that the light-like Wilson line in QCD without path 
ordering is given by 


Vexp[ig n dxf^A^^{x)T'^] 

Jxi 


1.X exp[-.9r‘[ 

dg 


2 I ■ D[A{xi)] 


I ■ 


d[gA{xi)]. 

dg 


(24) 
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which reproduces eq. ([5]) where the right hand side does not contain the path ordering V. 

Since the light-like Wilson line in QCD does not depend on the path but depends only 
on the end points jd , 5 ] we find from eqs. ( 125 ]) and (12 2 p that the non-abelian phase or the 
gauge link in QCD without path ordering is given by 


which reproduces eq. (Ej) which is used to study factorization of soft and collinear divergences 
in QCD where the right hand side of the above equation does not contain the path ordering 
V. 


III. NON-ABELIAN GAUGE TRANSFORMATION OF LIGHT-LIKE WILSON 
LINE IN QCD WITHOUT PATH ORDERING 


In order to study the gauge transformation of the light-like Wilson line in QCD without 
path ordering we proceed as follows. The non-abelian gauge transformation is given by 

TM “ n ) = U(x)T‘Al(x)U-\x) + -[dM^)]U-\x) ( 26 ) 

'^9 

where 


U{x) = 


( 27 ) 


Since the matrices are non-commuting we find from eq. 


T'^U-\x) = = r“[l + {-ig)T^u\x) + ^^T^T^u\x)u%x) 




2 ! 


+ 


-tgf 


-ig) 


T>’T'^T'^u;\x)u%x)u%x) + ^^I^T'^T‘^T'^T^uj\x)u;‘^{x)uj%x)u^{x) + ...].( 28 ) 


3! V 2 V 2 V 2 41 

By repeated use of the commutation relation 


j'ya ^ jabcrj^c 


(29) 


we find from eq. 


-ig? 


T^U-\x) = [T“ + {-ig)T’’u\x)T'^ + ^^^T^T^u;\x)uj%x)T^ 


3! 


4! 
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2 ! 


-tgf 


igf 


2 ! 


T^u\x)T^uj‘^{x)ir'^^u\x)T^ + ^^T^uj\x)ir‘^^u^{x)ir‘^^u\x)T<^ 


2 ! 


^^^\r'^^u\x)ir^<^u%x)iP'^^u<^{x)T’^ + ^^-^T’^uj\x)T'^uj%x)T'^u<^{x)ir^<^co^{x)T<^ 


3! 

-^gr 


3! 


^gY 


2 ! 2 ! 
-tgY . 


T^uYx)T^uj‘^{x)ir‘^^uj‘^{x)iP^^uj%x)T^ + ^^T’’uj\x)ir^^uj^{x)iP^^u\x)if^^u%x)T 


3! 


ir’’>^u\x)iP‘^^u'^{x)if'^^u‘^{x)ir^u^{x)T^ + ...] 


4! 


which gives after simplihcation 


T^U-\x) = [T“ + {-ig)T^uj\x)T'^ + 


-igf 


■^aTT-lt^\ _\ma , ^ , ..^rT.b. .b(^\rT^a , ^ -T'^T^Oj\x)uj\x)T° 


{-igf 

3! 


2 ! 

[-wf 


■T^rT‘^uj\x)uj%x)uj'^{x)T‘^ + ^-f^TV(a:)r"a;"(x)rV(x)T"a;"(x)T“ + ... 


4! 


-igf 


+ [1 + {-ig)T'^uj\x) + ^^^T'^u\x)T’^oj%x) + ..]{-ig)ir<^^u'^{x)T^ 


2 ! 

^igf 


^g?.. 


+ [1 + {-ig)T'^uj\x) + ^^;^T^u\x)T^oj\x) + ..]^^;^ir‘^oj^{x)if^^uj^{x)T 


2 ! 

{-^g? . 


2 ! 


+ [1 + {-ig)T^uj\x) + ..]^^^ir^^oj\x)ir3u\x)if^‘^^oj\x)T^ 


3! 


-igf . 


+[i + ..]^^ir^uj\x)ir'^u\x)if^^^uj\x)ir^u\x)T^ +...]. 


4! 

From eq. fl3T]) we hnd 


(31) 


T^U-\x) = U-\x)[T^ + {-ig)irV{x)r + 


-igf ■ 


2 ! 


if^9d^9^x)if’^^u’^{x)T^ 


-tgf . 


3! 

-igf ■ 


ir^^ufx)ir^>^uj%x)ip‘^'^u\x)T’^ 


4! 


ir^^uj\x)iP^’^u%x)if‘^^u\x)iriuj%x)T^ + ...] 


(32) 


which gives 

U{x)T^U-\x) = [e-<^^^%bT’’ (33) 

where Mab{x) is given by eq. ([8]). From eq. fl33|) we hnd 

f/(x)T“A“(x)f/-i(x) = (34) 

Similarly by simplifying inhnite nnmber of non-commnting terms in [dfj_U{x)]U~^{x) we 
hnd je] 


1 pgM{x) _ 1 

-\d,u{x)]u-\x) = |a,‘^‘(x)]r" 


(35) 


(30) 
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where Mah{x) is given by eq. ([8]). 

Hence by using eqs. fIMl) and fl3^ in eq. fl26|) we find 

pgM{x) _ 1 

A';(x) = (36) 

which is the finite gauge transformation in QCD where Mab{x) is given by eq. (IH]). Under 

infinitesimal gauge transformation we find from eq. (l36il 

+ gr’^^u%x)A^^\x) + d>^u^{x) (37) 

which is the infinitesimal gauge transformation in QCD which is familiar in the literature 

Q, 

When is the SU(3) pure gauge we find by using eq. ([7]) in fl36l) that 

2gM{x) _ 1 

A"“(x) = (38) 

By using eq. flTTll in fl38ll we find 

2g^N(x) _ 1 

U (39) 

where Nab{x) is given by eq. flT^ which is independent of g because /3“(a;) is independent 
of g, see eq. flTT]) . Multiplying the matrix gN[x) in eq. fl59]) we obtain 

D'^[y4'(x)]/3“(x) = [d^^\x)] (40) 

where 

D“^[W(x)] = + gr^A%{x). (41) 

By multiplying g in eq. fl3^ and then taking the derivative with respect to g we find 

= 49 \d“l3\x)]. (42) 

Using eq. fl40|l in fl42|) we obtain 

'^^AA^^igD^[A'(x)\li-‘(x). (43) 

By using eq. flTT]) in (I43|) we find 

d[gA'^^{x)] ^ ^ d^[A{x)]u\x). (44) 
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By multiplying the same independent four vector in eq. we obtain 


I ■ = 41- D[A{x)]u\x). 

cLg 

By dividing I ■ D[A'{x)] from left in eq. fH5|) we find 

■ ‘ ;. iMMi- = 2^“(.). 


( 45 ) 


2 I ■ D[A'{x)] dg 
Under the non-abelian gauge transformation as given by eq. (|26|) we find from eq. 

Hence from eqs. fH7|) . fH 6 |) . (l22i) and fl27|) we find 


(46) 


(47) 


$'(a;) = U{x) $(x), 


<h''^(x) = <h^(x) U ^(x) 


(48) 


which is the gauge transformation of the non-abelian phase in QCD under the non-abelian 
gauge transformation as given by eq. fl26|l . 

From eqs. fl22i) . fl2^ and (1481) we find 


d\i-A'‘^{x+i\)T‘^ = U{x)Ve~^^ , U{x) = ^ 49 ) 


which is the gauge transformation of the non-abelian gauge link in QCD under the non- 
abelian gauge transformation as given by eq. (|26ll . 

From eqs. (I23l) and fj4^ we find that, under the non-abelian gauge transformation as 
given by eq. fl26|l . the light-like Wilson line in QCD transforms as 


Ve 


igj^f dx^A'-{x)T- ^ 


Ve 


ig dxi^A'^{x)T° 


U ^(Xi), 


u{x) = ( 50 ) 


IV. EMERGENCE OF LIGHT-LIKE WILSON LINE IN QCD IN THE PROOF 
OF FACTORIZATION THEOREM AT HIGH ENERGY COLLIDERS 


Note that in we have shown that the light-like Wilson line in QCD naturally emerges 
when path integral formulation is used to prove NRQCD factorization at all order in coupling 
constant in heavy quarkonium production. Similarly, in we have shown that the light- 
like Wilson line in QCD naturally emerges when path integral formulation is used to prove 
factorization of soft and collinear divergences of the gluon distribution function at high 
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energy colliders at all order in coupling constant. In this paper we will prove that the 
light-like Wilson line in QCD naturally emerges when path integral formulation is used to 
prove factorization of soft and collinear divergences of the quark distribution function at 
high energy colliders at all order in coupling constant. 


The generating functional in the path integral method of QCD is given by 

Z[J,V,v] = 


16 


,ll7| 


e-/- 






where is the external source for the quantum gluon held Q^°‘{x) and 77 * (x) is the 

external source for the Dirac held of the quark and 


F‘[Q] = - dMx) + gf^CfMQKx), 


pa2 


flU 


[Q] = F»“'“|Q|F;„|Q|. (52) 


The light-like quark traveling with light-like four-velocity produces SU(3) pure gauge 
potential at all the time-space position x^ except at the position x perpendicular 

to the direction of motion of the quark (/ • x = 0) at the time of closest approach 2, 00. 
Hence the soft and collinear behavior of the non-perturbative correlation function in QCD 
due to the presence of light-like Wilson line in QCD can be studied by using path integral 
formulation of the background held method of QCD in the presence of SU(3) pure gauge 
background held 0 . 0 . 


Background held method of QCD was o rig inally formulated by ’t Hooft 18| and later 


extended by Klueberg-Stern and Zuber [l9|, [2^ and by Abbott [l7|. This is an elegant for¬ 


malism which can be useful to construct gauge invariant non-perturbative green’s functions 


in QCD. This formalism is also useful to study 


chromo held 2]| via Schwinger mechanism 22l. to compute (3 function in QCD 


perform calculations in lattice gauge theories 
constant in the presence of chromoheld 251. 


quark and gluon production from c 


assical 


23[, to 


24| and to study evolution of QCD coupling 


It can be mentioned here that in soft collinear ehective theory (SCET) 26| it is also 
necessary to use the idea of background helds [l7| to give well dehned meaning to several 
distinct gluon helds [9|. 

Note that a massive color source traveling at speed much less than speed of light can 
not produce SU(3) pure gauge held 0 0 7]. Hence when one replaces light-like Wilson 
line with massive Wilson line one expects the factorization of soft/infrared divergences to 


(51) 
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break down. This is in conformation with the hnding in m which used the diagrammatic 
method of QCD. In case of massive Wilson line in QCD the color transfer occurs and the 
factorization breaks down. Note that in case of massive Wilson line there is no collinear 
divergences. 


The generating func 
of QCD is given by [17 


ional in the path integral formulation of the background held method 


19| 


Z[A, J,ri,ri\ = J[dQ][d'ifj][d'ilj] det{ 


SG^iQ). 

5uj^ ' 




(53) 


where the gauge hxing term is given by 

G\Q) = = D,[A]Q^<^ 

which depends on the background held and 

+Q\ = d,[Ai + oa - +q ;]++ oy k+« i- 


(54) 


(56) 


We have followed the notations of 


17 


19| and accordingly we have denoted the quantum 


gluon held by <5^“ and the background held by 

Note that the gauge hxing term in eq. fl55]) [where G°‘{Q) is given by eq. 

(jHTll ] is invariant for gauge transformation of A)): 


= g+ 9^0;“, (type I transformation) 


provided one also performs a homogeneous transformation of Q® 



SQl = af'^Qy- 


(56) 


(57) 


The gauge transformation of background held ^4“ as given by eq. (l56|l along with the 
homogeneous transformation of <5“ in eq. fl571) gives 

6{A; + Q“) = gr’^^iAl + Ql)u^ + (58) 

which leaves + Q] invariant in eq. fl53|l . 

For hxed v4“, he., for 


= 0, (type II transformation) 


(59) 
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( 60 ) 


the gauge transformation of 





gives eq. fl58l) which leaves + Q] invariant in eq 

It is useful to remember that, unlike QED js], Ending an exact relation between the 
generating functional Z[J, g, fj\ in QCD in eq. (l5T|l and the generating functional Z[A, J, rj, g] 
in the background field method of QCD in eq. fl5^ in the presence of SU(3) pure gauge 
background field is not easy. The main difficulty is due to the gauge fixing terms which are 
different in both the cases. While the Lorentz (covariant) gauge fixing term —in 
eq. fl5T]) in QCD is independent of the background field the background field gauge 

fixing term —in eq. in the background field method of QCD depends on 
the background field y4^“(a;) where G'^{Q) is given by eq. flHTl) 17H19|. Hence in order to 
study non-perturbative correlation function in the background field method of QCD in the 
presence of SU(3) pure gauge background field we proceed as follows. 

By changing Q —)■ Q — H in eq. fl53D we find 


(61) 


Z[A,J,r],f]] = 


det(- 




where the gauge fixing term from eq. fIST)) becomes 


G}(Q) = = D^[A]Q^^‘^ - d^A^^\ (62) 

and eq. fl57)) [by using eq. fl56|l . type I transformation [^, 19|] becomes 


5Ql = gr^^Qlu^ + d,u\ (63) 

The eqs. (I62|) and fl63|) can also be derived by using type II transformation which can be 
seen as follows. By changing Q —)■ Q — H in eq. flS5]l we find eq. fIM]) where the gauge fixing 
term from eq. flST)) becomes eq. (16^ and eq. (16U]) [by using eq. flS^ ] becomes eq. (165]) . 
Hence we obtain eqs. (EH), dSZD and (l63|) whether we use the type I transformation or type 
H transformation. Hence we find that we will obtain the same eq. fl8T|) whether we use the 
type I transformation or type H transformation. 

The equation 


= QU^) + gr’^'^u;'^{x)Qlix) + d.co^ix) 


(64) 
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in eq. (163|) is valid for infinitesimal transformation (a; << 1) which is obtained from the 
finite equation 


= f/(x)T“Q“(x)17-'(x) + -[d,U{x)]U-\x), 


U{x) = (65) 


Simplifying infinite numbers of non-commuting terms in eq. fl6^ [by using eq. (jSS]) and 
Q] we find that 

, , pQ^ix) _ 1 

Q'lix) = U Mab{x) = r^^u%x). (66) 

Changing the variables of integration from unprimed to primed variables in eq. fl6Tll we 


Z[A, J, ri, fi] = [ [dQ'j [#'][#'] 


find 


Su^ 

This is because a change of variables from unprimed to primed variables does not change 
the value of the integration. 

Under the hnite transformation, using eq. fl6^ . we find 

ya 


m = [dQ] detl^l = |d0] det[[e»"l">]] = |d<J]exp|Tr(ln[e»“<*)])] = \dQ] (68) 


where we have used (for any matrix H) 


detif = exp[Tr(lnif)]. 


(69) 


The fermion field transforms as 


^'(x) = (70) 

Using eqs. (166]) and (176]) we find 

[#'][#'] = [#][#], -m + - m + 

= F^iAQi ( 71 ) 

Using eqs. (1681) and (1711) in eq. (167]) we find 
ZIA, J, p, g] = J [dQ] [#] [#] 

det(^^5i^) gi f rfU[-iF“2jQ]_^(G“(Q'))2+j.Q'+p[i7Ma^_m+gT“77‘Q“]y+i?y'+y'r)]_ ^^^2) 
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From eq. fl6^ we find 


GJiQ') = + gr'^AlQ''^ - a„/l'‘“. 


( 73 ) 


By using eqs. (1661) and ([7j) in eq. d73|) we find 

pgM(x) _ -| 

g;( 0 ') = Mm 

pgM(x) _ -| 

WMcM) + I ].4a.^‘‘(^)ll 


which gives 


-a„|a^u‘(i 


^gM{x) _ 

gM{x) 

^gM(x) _ ^ 


gM{x) 


_ ah 


( 74 ) 


G5(Q') = a"[|e»"W]„tQ‘(i)| 

r p9M(x) _ 1 


pgM(x) _ 1 

]iMUQi(x) + u [Mm- ( 78 ) 


be 


gM{x) 


From eq. (j75ll we hnd 

GKQ') = ff'iim'rni.QUm + gr'"M(x) 

which gives 


^gM(x) _ ^ 

gM{x) 


]lM'UQi(x)] ( 76 ) 


be 


G}{Q') = [e^^^%bd^Ql{x) 

+Q’l{x)d>^[[e^^^%b] + [d>^u^{x) 


^gM{x) _ ^ 

gM{x) 


wmimmMx)]- (77) 


- be 


From 


6| we find 


= ig[d^u\x)] 


^gM(x) _ 

gM{x) 


M) = Mix) 


ab 


( 78 ) 


which in the adjoint representation of SU(3) gives (by using Tj“ = 

' ^gM{x) _ 


= [d>^u%x)] 

Using eq. (17^ in d77|) we find 


gM{x) 


gf^c^e^i^,, M,,(x) = r’^'^co'^ix). (79) 


be 


GUQ') = M’UMJx) 


(80) 
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which gives 


(Gj(Q')y = (a„<?'‘“(x))‘ 


(81) 


Since for n x n matrices A and B we have 


det(x4i?) = (dety4)(det5) 


(82) 


we find by using eq. flHOD that 


= det|-] = det[|e» *]., J . 


det[[e' 


,gM(x)l 


det[ 


5{d^^Q’l{x)) 
5u^ ^ 


= exp[Tr(ln[e5^(*)])] det[ 


Su^ 


= det[ 


d{df,Qf^^{x)). 

du’’ 


(83) 


Using eqs. fl8T|) and fl83|) in eq. fl72D we find 
Z[A,J,r],f]] = 




From eqs. ([7]) and fIBB]) we find 


Q"(x) - A-M = |e»"'“’>|.xQt(x), 




(85) 


Note that eqs. flM)l and fl57]l are valid whether we use type I transformation [eqs. flSB]) and 
flSTjl ] or type II transformation [eqs. fl59|l and fl60|l ]. 

Since we have used eq. (I26|) to study the gauge transformation of the Wilson line in QCD 
we will use type I transformation, see eqs. fl56|) and fl57)) . in the rest of the paper which 
gives for finite transformation jw. 19[ 


= [e^^^%Mx), 


MM) = 


( 86 ) 


From eqs. 


and fl86l) we find 


Z[A,j',rj,rj] = 
5u^ 
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Under the non-abelian gauge transformation the fermion source transforms as 17|, Il9 | 

( 88 ) 

From eqs. fl70|l and fl88|l we hnd 


-ijj'r]' = tfjri 


(89) 


which gives from eq. flHT)) 

Z[A,J',r]',f]'] = J [dQ][d'i(j][d'ijj] 

5uA 

Hence from eqs. and fl5Tl) we hnd that in QCD 

Z[J,r^,rf^ = Z[A,J\i,f^'] (91) 

when the background held H^“(x) is the SU(3) pure gauge as given by eq. ([3]). The corre¬ 
sponding relation in QED is given by 




(92) 


when the background held A^{x) is the U(l) pure gauge as given by eq. ([T]). Note that, 
unlike eq. fl9T]l in QCD there is no J' in eq. fl92l) in QED because while the (quantum) 
gluon directly interacts with classical chromo-electromagnetic held the (quantum) photon 
does not directly interact with classical electromagnetic held. 

die non-perturbative correlation function of the type < O|'0(a;)'0(x')|O > in QCD is given 

8 | 


by 


< Q\'il){x)'ip{x')\Q >= 


-7T^[<^,h,h]|j=r,=f)=0- 


(93) 


6ri{x) 5rj{x' 

Similarly the non-perturbative correlation function of the type < O|'0(a;)'0(a;')|O in the 
background held method of QCD is given by 


< 0\i>{x)ij{x')\0 >A= -^^-^^Z[A,J,r],fi\j=^=f^=o- (94) 

When background held H^“(a;) is the SU(3) pure gauge as given by eq. (|3]) we hnd from eqs. 
(1911). (193]), (EID, (IH6D and (JEEj) that 

< O|'i^(x)'0(a;')|O >=< 0|'i^(x)<F(a;)<F^(x')'^(x')|0 >a (95) 
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which proves factorization of soft and collinear divergences at all order in coupling constant 
in QCD where [see eq. fl2^ and jd, 5|] 

(96) 


roo 

d)(x) = Pexp[-z^T“ / dXl ■ A^{x + /A)] = 

Jo 


is the non-abelian phase or the gauge-link in QCD. 

From eq. fl9^ we find that the correct definition of the quark distribution function at high 
energy colliders which is consistent with the number operator interpretation of the quark 
and is gauge invariant and is consistent with the factorization theorem in QCD is given by 


fq/pi^) = ^Jdy e X < P\tl>(0,y ,OTh*lPexpligT‘dz A*'‘(0,z , 0T)|]iii(0)|P > 

(97) 


gauges, in general 
. In eq. (1771) the 


which is valid in covariant gauge, in light-cone gauge, in general axia 
non-covariant gauges and in general Coulomb gauge etc. respectively 

is the Dirac field of the quark and A^°-{x) is the SU(3) pure gauge background field as 
given by eq. ([3]). 

Hence we find from eq. (1971) and from 4, Isj that the light-like Wilson line in QCD 
naturally emerges when path integral formulation of QCD is used to prove factorization of 
soft and collinear divergences at all order in coupling constant in QCD processes at high 
energy colliders. 


V. CONCLUSIONS 

Unlike the Wilson line in QED the Wilson line in QCD contains path ordering. In this 
paper we get rid of the path ordering in the light-like Wilson line in QCD by simplifying 
all the infinite number of non-commuting terms in the SU(3) pure gauge. We have proved 
that the light-like Wilson line in QCD naturally emerges when path integral formulation of 
QCD is used to prove factorization of soft and collinear divergences at all order in coupling 
constant in QCD processes at high energy colliders. 
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